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. i-H . 1 INTRODUCTION 



ABSTRACT 

The early X-ray light curve of Gamma Ray Bursts (GRBs) is complex, and shows a typical 
steep-flat-steep behaviour. The time T a at which the flat (plateau) part ends may bear some 
important physical information, especially if it plays the same role of the so called jet break 
time tj Ct . To this aim, stimulated by the recent analysis of Willingale et al., we have assembled 
a sample of GRBs of known redshifts, spectral parameters of the prompt emission, and T a . 
By using T a as a jet angle indicator, and then estimating the collimation corrected prompt 
energetics, we find a correlation between the latter quantity and the peak energy of the prompt 
emission. However, this correlation has a large dispersion, similar to the dispersion of the 
Amati correlation and it is not parallel to the Ghirlanda correlation. Furthermore, we show 
that the correlation itself results mainly from the dependence of the jet opening angle on the 
isotropic prompt energy, with the time T a playing no role, contrary to what we find for the 
jet break time ij Ct . We also find that for the bursts in our sample T a weakly correlates with 
-Ey.iso of the prompt emission, but that this correlation disappears when considering all bursts 
of known redshift and T a . There is no correlation between T a and the isotropic energy of the 
plateau phase. 

Key words: gamma rays: bursts, X-rays: general, radiation mechanisms: general. 



X 



Before Swift, the afterglow observations were typically starting 
only several hours after the burst, when the flux typically showed 
a single power law decay. The detection of GRB afterglows pro- 
vided strong confirmation for the fireball shock model (Rees & 
Meszaros 1992), which explains it in terms of forward shock emis- 
sion running into the circumburst medium. In many of the well 
sampled afterglows, the steepening in the light curve at late times 
(~ 10 5 s) has been attributed to a narrow conical jet of semiaper- 
ture opening angle 0j whose edges become visible as it decelerates 
and widens (e.g. Rhoads et al. 1997). The deceleration is due to the 
expansion of the GRB outflow into the interstellar medium (ISM) 
whose properties are still poorly understood. By assuming a uni- 
form density ISM the widely dispersed isotropic energy -E 7 ,i so of 
GRBs considerably clusters if corrected for the collimation factor 
/ = (1 - cosflj) (Frail et al. 2001). Later, Ghirlanda, Ghisellini 
& Lazzati (2004; hereafter GGL04) discovered that the collimation 
corrected GRB energies _E 7 = E lt i so f are tightly correlated with 
the rest frame peak energy E p of the vF v prompt spectrum (so 
called "Ghirlanda" correlation). Nava et al. (2006; hereafter N06) 
discovered that, in the case of a wind-like circumburst density pro- 
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file (which should be more likely if the progenitor is a massive 
star), this correlation still exists and becomes linear. Liang & Zhang 
(2005) discovered the phenomenological counterpart of these cor- 
relations showing that a tight correlation exists between the -E 7 ,i so , 
E p and tj ct (observed in the optical light curves) computed in the 
source frame. N06, indeed, demonstrated that this phenomenologi- 
cal correlation is consistent with the E p -E-, correlations computed 
either in the case of a homogeneous and wind medium (HM and 
WM, respectively). 

The jetted-fireball model predicts that when the bulk Lorentz 
factor F ~ due to the deceleration of the expanding fire- 

ball, an achromatic break appears in the afterglow light curve. Such 
achromaticity is what distinguish a jet break (i.e. due to the geom- 
etry of the event) from a spectral break due to the time-evolution 
of the characteristic frequencies of the spectrum. Most of the jet 
breaks of the pre—Swift GRBs were, indeed, achromatic, but in a 
narrow frequency band, from the near-IR to the optical. Moreover, 
the typical jet breaks were observed between several hours to few 
days after the GRB trigger: in order to identify this transition in the 
optical light curves a systematic follow up of the optical afterglow, 
up to several days after the jet break, was required. 

It was highly expected that the multi-wavelength Swift satel- 
lite (Gehrels et al. 2004) would clearly detect achromatic X-ray to 
optical breaks. Instead, Swift discovered that the X-ray light curve 
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of GRBs is much more complex than thought (e.g. Burrows et al. 

2005) . The "canonical" X-ray light curve disclosed by the Swift ob- 
servations is composed by an initial steep decay phase (commonly 
interpreted as the prompt emission tail due to off axis radiation of 
a switching-off fireball) followed by a shallow decay phase. At the 
end of this second phase there is a break that marks the beginning of 
the "normal" decay phase, observed also before Swift (i.e. when the 
X-ray follow up - mostly with BeppoSAX - started, several hours 
after the trigger). In addition, Swift found that the X-ray afterglow 
is characterised by several early and late times flares superposed to 
the typical power law flux decay (Nousek et al. 2006, O'Brien et al. 

2006) . In several Swift bursts the optical light curve does not track 
the X-ray one while in some cases flares and multiple breaks have 
been observed in the optical too. Due to the richness and complex- 
ity of the light curves, the identification of a break as a jet break 
requires care (Ghirlanda et al. 2007 - G07 hereafter). 

Some interpretations of the complex multi-break X-ray light 
curve of Swift GRBs have been proposed: while it is commonly 
accepted that the early steep decay can be interpreted as due to 
the curvature effect of a switching-off fireball (see e.g. Kumar & 
Panaitescu 2000), the origin of the intermediate shallow decay is 
still controversial. It has been proposed (Ghisellini et al. 2007) that 
flat X-ray could be "late prompt" emission due to a late central 
engine activity producing shells with a time-decaying bulk Lorentz 
factor r. In this scenario the flat-to-steep transition occurs when T 
of these late shells becomes ~ l/8j . 

Recently, Willingale et al. (2007; hereafter W07) interpreted 
the steep-flat-steep X-ray light curves as due to the superposition 
of two separate and independent components: the prompt emission 
(interpreted as internal dissipation) produces the first steep decay 
while the afterglow emission (external dissipation) is responsible 
for the shallow and the final steep components. The modelling, in- 
troduced by W07, identifies an early break (typically ~ 10 3 -10 4 
s) in the X-ray light curve, called T a , that marks the shallow-to- 
steep transition. W07 show that if T a is used (with the same for- 
malism adopted for tj Gt in the HM case) to compute an angle 0j, a , 
the collimation corrected energy -E 7 , a = -E 7 ,i so (l — cos#j, a ) is 
correlated with E p . W07 argue that this correlation has the same 
slope of the E p — E 1 (Ghirlanda) correlation (found with the jet 
break time ij c t). The two correlations are simply displaced by a 
factor ~ 26 in _E 7 . This implies that the two break times T a and 
ijet are proportional (i.e. tj ot ~ 90T a ). In the sample of W07 only 
GRB 050820A has both T a and tj e t measured, and indeed in this 
case tj e t~ 100T a . W07 conclude that if there is a link between T a 
and tj c t, it is doubtful that both times are actually jet break times. 
However, W07 stress that it is possible that the end of the plateau 
phase (corresponding to T a ) depends on the total energy of the out- 
flow. They also point out that it would be interesting to study the 
phenomenological counterpart of the E p — i? 7 . a correlation simi- 
larly to what has been done with the Liang-Zhang correlation with 
respect to the Ghirlanda one. 

The scope of this paper is to investigate the role of the break 
time T a in defining the E p — i? 7 , a correlation and if a Liang- 
Zhang like correlation exists between the quantities E p , E lt i so and 
T a (where E p and T a are expressed in the source rest frame). To 
this aim, we enlarge the sample used by W07 from 14 to 23 GRBs 
and derive the collimated corrected energy S 7 , a using T a as a jet 
break and we reproduce the correlation between E p and i? 7 , a and 
test, through the scatter of these correlations, the nature of T a . 

We adopt a standard cosmology with fl m = 0.3, S1a = h = 

0.7. 



2 THE SAMPLE 

In order to build the E p — _E 7 , a correlation we need the redshift z, 
the peak energy of the prompt emission spectrum E p and, from 
the X-ray light curve, the time T a . We found, up to November 
2006, 16 objects with all these quantities published in the literature. 
We added 7 more GRBs (050126 050505 050908 060223A 060522 
060607A 060908) for which we analysed the Swift-BKT (15-150 
keV) spectrum and could satisfactorily constrain the spectral peak 
energy E p fitting the spectrum with a cut-off power law model. In 
these cases we also tested our results (by the _F-test) finding that 
the fit with a cut-off power law is better than a single power law. 
The details of the fitting, together with the results for other Swift 
bursts, will be given elsewhere (Cabrera et al. 2007 in preparation). 

The sample of W07 which is used to build the E p — i? 7 , a 
correlation contains 14 long GRBs. However, it is hard to identify 
those bursts present in both samples because W07 do not give the 
names of those bursts they use to compute the E p — -E 7 , a correla- 
tion. 

Our sample is composed by 23 GRBs whose properties are 
listed in Table Q] In most cases the value of T a is taken from the 
compilation of W07, except for two cases: 

• for GRB 050922C we have found that the X-ray light curve 
does not clearly show the shallow part and, therefore, the value of 
T a is hardly constrained. A still acceptable fit is obtained by setting 
T a ~10 4 seconds (a factor 26 larger than what derived by W07 - 
see G07); 

• for GRB 060124 the fit of the X-ray light curve should ac- 
count for the late time achromatic break (Romano et al. 2006; Cur- 
ran et al. 2006) which corresponds to the jet break time (G07). With 
the inclusion of this further late time break T a ~ 5000 s (a factor 8 
smaller than what reported in W07). The fit is shown in Fig.Q] 

Moreover, since the sample of W07 comprises GRBs up to 
August 2006, for GRB 060908, GRB 060927 and GRB 061121, 
we have estimated T a following the procedure outlined in W07, 
and using the publicly available X-ray light-curve^]. The corre- 
sponding light curves and fitting models are shown in Fig.Q] 

In the following sections we analyse the E p — £ 7 , a correla- 
tion (found using T a as an angle estimator) and compare it with the 
E p — E-y correlation (found using tj ct as angle indicator). We also 
study the relevance of either T a and tj ct in "collapsing" the scatter 
of the Amati correlation (Amati et al., 2002; Amati 2006), in the 
respective collimation corrected correlations. 

For the E p — i? 7 , a correlation we use the sample of 23 GRBs 
of the Swift-em reported in Tab. 1. For the comparison with the 
E p — i? 7 correlation we use the most updated sample of 23 bursts 
for which ij ot has been firmly identified from the optical light curve 
reported in Gofl Within this sample (with ij ot ) there are 6 GRBs 
of the Swift-em for which we also have an estimate of T a . To re- 
produce the Amati correlation (which only requires the knowledge 
of the redshift and of the spectral parameters) we have used those 
GRBs that appear at least in one of the two samples. This larger 
sample contains 40 GRBs. 



1 http://astro.berkeley.edu/~nat/swift/ by N. Butler, see also Butler & Ko- 
cevski (2007) 

2 In the present analysis we do not include the lower limits on tj et dis- 
cussed in G07. 
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GRB 


z 


E p 




Instr" 


Ref 6 




Ref d 








keV 


erg 






s 






050318 


1.44 


115±27 


2.03±0.34 


SWI 


1 


2.01 [1.37-2.68] 


W07 




050401 


2.90 


501±117 


41±8 


KON 


2 e 


3.87 [3.69-4.11] 


W07 




0504 16A 


0.653 


28.6±8.3 


0.083±0.029 


SWI 


3 


3.19 [2.69-3.48] 


W07 




050525A 


0.606 


127±5.5 


2.96±0.64 


SWI 


4 


2.92 [1.69-3.06] 


W07 




050603 


2.821 


1333±107 


59.8±4.0 


KON 


5 


4.83 [3.69-5.13] 


W07 




050820A 


2.612 


1325±277 


97.8±8.0 


KON 


6 


3.96 [3.84^1.04] 


W07 




050922C 


2.198 


417±118 


4.61±0.87 


HET 


7 


2.58 [2.48-2.69] 


W07 
















3.98 [3.68^1.16] 


G07 




051109A 


2.346 


539±381 


7.58±0.94 


KON 


8 


3.93 [3.70^1.08] 


W07 




060115 


3.53 


281±82 


8.0±1.5 


SWI 


9 


3.86 [2.86-5.47] 


W07 




060124 


2.297 


636±162 


43±4.0 


KON+SWI 


10 


4.60 [4.47^1.69] 


W07 
















3.70 [3.40-3.88] 


this paper, 


see Fig.[T] 


060206 


4.048 


381±98 


4.75±0.76 


SWI 


11 


3.86 [3.68^1.00] 


W07 




060418 


1.489 


572±114 


12.8±1.1 


KON 


12 


3.44 [3.26-3.58] 


W07 




060526 


3.21 


105±21 


2.58±0.26 


SWI 


13 


3.84 [3.51^4.31] 


W07 




060707 


3.43 


292±71 


6.83±1.49 


SWI 


14 


3.58 [2.87^1.05] 


W07 




060927 


5.6 


473±116 


9.69±1.62 


SWI 


15 


3.60 [3.30-3.78] 


this paper, 


see Fig.[T] 


061121 


1.314 


1289±153 


26.1±3.0 


KON/RHE 


16-17 


3.90 [3.60-4.08] 


this paper, 


see Fig.[T] 


050126 


1.29 


202±49 


1.10±0.30 


SWI 


this paper 


2.34 [1.34-5.64] 


W07 




050505 


4.27 


622±211 


19.5±3.1 


SWI 


this paper 


4.39 [4.15^1.87] 


W07 




050908 


3.344 


191±40 


2.05±0.34 


SWI 


this paper 


3.31 [2.03^1.17] 


W07 




060223A 


4.41 


341 ±69 


4.54±0.72 


SWI 


this paper 


2.73 [2.37-3.05] 


W07 




060522 


5.11 


415±78 


8.02±1.65 


SWI 


this paper 


2.86 [2.61-3.32] 


W07 




060607A 


3.082 


535±164 


12.2±1.8 


SWI 


this paper 


4.75 [4.73^1.78] 


W07 




060908 


2.43 


487±116 


9.17±1.57 


SWI 


this paper 


3.10 [2.80-3.28] 


this paper, 


see Fig.[T] 



Table 1. GRB sample used to study the E p — S 7 , a correlation. The peak energy E p is in the source rest frame. The time T a is in the observer frame. i? 7 ,i so 
is given in units of 10 52 erg. "The + sign indicates that the spectral parameters are derived from the joint spectral fit of two different instruments. A sign 
/ indicates that the reported values of B p and i? 7 ,i so are an average of the values found by two different instruments. The instruments are: SWI=Swfyf/BAT, 
KON=Konus/Wind, HET=HETE/Fregate, RHE=RHESSI. ''References for the spectral parameters: (1) Perri et al., 2005; (2) Golenetskii et al., 2005a; (3) Sato 
et al., 2006; (4) Blustin et al., 2006; (5) Golenetskii et al., 2005b; (6) Cenko et al., 2006; (7) Crew et al., 2005; (8) Golenetskii et al., 2005c; (9) Barbier et al., 
2006; (10) Romano et al., 2006; (11) Palmer et al., 2006; (12) Golenetskii et al., 2006a; (13) Schaefer 2007; (14) Stamatikos et al., 2006a; (15) Stamatikos et 
al., 2006b; (16) Golenetskii et al., 2006b; (17) Bellm et al., 2006. c Values in square brackets are the lower and the upper 90% confidence limits. ^References 
for T a : Willingale et al. 2007 (W07); Ghirlanda et al. 2007 (G07). e The values of E p and i? 7 ,i so derive from the average of the two spectral fit presented in 
the corresponding reference. For GRB 050922C and 060124 we report the value of T a given in W07 (first line) and the value obtained in G07 and this paper 
(second line), respectively. 

where cLl is the luminosity distance, is the 7-ray fiuence in the 
observed energy band, k is the bolometric correction factor needed 
to find the energy emitted in a fixed energy range (here, 1-10 4 
keV) in the rest frame of the source. For most pre-Swift bursts the 
spectrum is fitted with the Band function (Band et al. 1993) com- 
posed by two smoothly joint power laws. On the other hand, several 
Swift bursts, due to the narrow energy range of the BAT instrument 
(15-150 keV), have a spectrum best fitted with a cutoff-power law 
model. As discussed in Firmani et al. (2006), in these cases the 
extrapolation of the cutoff-power law model up to 10 4 keV may 
considerably underestimate the energy, if the spectrum is, instead, 
a Band function. To account for this limitation of the spectral fits 
in the Swift era, the value of i? 7 ,i so reported in Tab. 1 is an aver- 
age between the value derived with the cutoff-power law model 
and that derived assuming a Band model with high energy photon 
index j3 = —2.3 (see also G07). 



3.1 Correlation analysis 

For the 23 GRBs of Tab. 1 we computed the angle #j, a , where the 
subscript "a" means that the angle has been estimated using T a as 
a jet break time, and the collimation corrected energy -E 7 , a in both 
the WM and HM case. The density in the homogeneous case is 



3 ANALYSIS AND RESULTS 

First we study the Ghirlanda like correlation and the Liang & Zhang 
like correlation defined using T a in place of tj ct . We use the formal- 
ism defined in N06 where, for a generic time break ibr Ca k, the jet 
opening angle is defined afl: 




where n represents the circumburst density in the homogeneous 
case (HM) and A* is the normalisation of the wind (WM) density 
profile (n(r) oc A-,.r~' 2 - see N06 for more details). Here ?7 7 is the 
radiative efficiency of the prompt phase assumed to be 0.2 for all 
bursts. 

In both cases of Eq. [TJ it is necessary to know the isotropic 
bolometric rest frame equivalent energy i5 7 ,i so . This is 

1 + Z 

3 We use the notation Q x = lO^Q, using cgs units. 
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Figure 1. The early BAT and late XRT light curves of GRB 060124, GRB 060908, GRB 060927 and GRB 061 121 modelled with the function proposed by 
W07. The two model components are shown by the long-dashed lines while their sum is shown by the solid line. The fits were performed excluding the flares 
of the light curves. The long dashed vertical lines show the T a in each light curve (for GRB 060124 we also report the jet break time found in the optical light 
curve - see text). 



known only for some GRBs and for the others we take the value 
n = 3 ± 2.76; the density in the wind case is assumed equal for all 
GRBs and without error (see N06 and G07 for details). 

We report also the correlation defined with the jet break time 
tjet and recently updated with the Swift bursts in G07. Here we 
use the sample of G07 (Tab. 1 in that paper) excluding bursts with 
only a lower limit on tj ct . For comparison we also recomputed the 
E p — -E 7 ,i so (Amati) correlation with the sample of 40 GRBs. 

The fit is performed using the routine fitexy (Press et al. 1999) 
which weights for the errors on the involved variables. All the fit 
results can be expressed in the synthetic form: 

where -E 7 , a must be replaced by _E 7 or -E 7 ,i so according to the cor- 
relation considered. Here -Ebar represents the barycentre of the en- 
ergy values of the sample considered. 

The three correlations are shown in Fig.[2](homogeneous case) 
and in Fig.[3](wind case), and the fit results are listed in Tab.[2] For 
every correlation we report also the reduced \ 2 and the la scatter 
of the data points around the best fit line, modelled with a gaussian 
distribution. 

From the comparison of the results reported in Tab. 2 and in 
Figs. [2] and [3] we note: 



(i) the correlation E p — E ll& derived using T a as a jet break 
time is not parallel to the E p — _E 7 correlation (derived using 
tjet)- In particular, independently from the specific density pro- 
file considered, the E p — -E 7 , a correlation is much flatter than 
the E p — E-, correlation and in the homogeneous case it is even 
slightly flatter than the Amati correlation. By reconstructing the 
sample of W07 (through the quoted GCNs and the values given in 
that paper), and analyzing this smaller sample, we indeed obtain a 
E p — -E 7 . a correlation parallel to the E p — _E 7 one, although not 
as tight. It is the inclusion of the additional GRBs that changes 
the slope of E p — -E 7 , a correlation and further increases its scatter. 
Also the 6 GRBs with both T a and tj ot (encircled symbols) seem 
to define a E p — -E 7 , a parallel to the E p — E-, correlation, and this 
issue will be discussed in Sec. 3.2. 

(ii) the fit of the E p — -E 7 , a correlation (in terms of x 2 ) is signif- 
icantly worse than that of the E p — E-, correlation, which, instead, 
is statistically acceptable; 

(iii) the E p — -E 7 , a correlation has a scatter which is a factor 2 
larger than that of the E p — _E 7 correlation and is similar to that of 
the Amati correlation. 

Item (i) suggests that the relation between ij e t and T a , if any, 
is not trivial (i.e. tj e t is not simply proportional to T a as suggested 
by W07). The large scatter of the E p — -E 7 , a correlation, which is 
comparable to that of the E p — -E 7 ,i so correlation, raises the ques- 
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10 48 10 50 1 Q 52 1 54 

Energy [erg] 

Figure 2. Homogeneous circumburst medium. From right to left we show the E p — i? 7i i so correlation (filled squares - 40 GRBs), the E p — E-y correlation 
(filled circles - 23 bursts from G07) and the E p — B-y, a correlation (filled stars - 23 GRBs of Tab. 1). The shaded regions (for the three correlations) represent 
the 1, 2, 3<r scatter of the data points (computed perpendicular) with respect to the best fitted correlation (solid lines). Circled points (on the three correlations) 
represent the 6 Swift bursts for which both T a and tj ct are known. 




10 48 1 Q 50 1Q 52 1Q 54 

Energy [erg] 

Figure 3. Wind circumburst medium. Symbols as in Fig|2] 
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Correlation Density N K s i?bar X? Scatter 

profile (ct k ) (<r s ) (erg) 



Bp 


— f 7,a 


homog. 


23 


4.71 (0.37) 


0.55 (0.05) 


3.11 x 10 49 


1.95 (21 d.o.f.) 


0.16 




— ^7,a 


wind 


23 


4.95 (0.38) 


0.78 (0.07) 


4.14 x 10 49 


1.95 (21 d.o.f.) 


0.16 




— .E 7 


homog 


23 


2.89 (0.16) 


0.71 (0.04) 


4.14 x 10 50 


1.14(21 d.o.f.) 


0.09 


Bp 


— .E/ 7 


wind 


23 


3.21 (0.17) 


1.06(0.07) 


2.47 x 10 50 


0.91 (21 d.o.f.) 


0.08 


£ P 


— R ■ 




40 


3.77 (0.11) 


0.62 (0.02) 


1.31 x 10 53 


5.70 (38 d.o.f.) 


0.16 


£ P 


— R ■ 




23 


4.83 (0.21) 


0.65 (0.03) 


1.22 x 10 53 


2.89 (21 d.o.f.) 


0.15 



Table 2. Best fit results of the correlations involving the rest frame peak energy and the collimation corrected energy computed either using T a as a jet break 
time and the real jet break time tj ot . We also report the Amati correlation for the considered 40 GRBs and for the 23 defining also the E p — £7 7 ,a correlation. 
The collimation corrected correlations are reported for both the homogeneous and wind medium scenario. K and s represent the normalisation and slope, 
respectively, of the best fit correlation (see Eq.|5J. The lc scatter of the data points around the best fit correlation is also reported. 



tion whether the new added variable T a can be treated as a jet break 
time or not. If not, we are still left with the possibility that a com- 
pletely phenomenological correlation between E p and i? 7 ,i so and 
T a exists (where T a is the time in the source rest frame). This pos- 
sibility will be explored in the next section. 



3.2 The E p -E ltiso -Ti correlation 

In this section we study the E p -E^ t i SO — correlation. First we fit 
the multi-variable correlation with the 23 GRBs of Tab. 1 through 
a least-square method without considering the errors on the three 
variables. The result is: 



£ 7 ,i ao ,52 = (7.29 ± 0.; 



T 



0.012 s 



9) (s57kev) 

0.06±0.59 



(4) 




Note that the exponent of T a is very close to zero and this implies 
that T a does not reduce the scatter of the E p — E-y^o correlation. 
Indeed, the above result is almost identical to the E p — -E 7 ,i S o cor- 
relation defined by the 23 GRBs reported in Tab. 1 (see also the last 
row of Tab. 2). 

For consistency with the fits reported in the previous sections 

we also fitted the E p -E~. i so -T, correlation by weighting for the jcjuriisin.. u- . % /-> 1 *u „ c .u 

v '' a job defined by bq. |5J (bottom histogram). One can see that the scatter of the 

errors on the three variables. In this case we find a different result: 



Figure 4. Scatter distribution of the 23 GRBs of Tab. 1 around the best 
fit E p — E-v i so correlation (top histogram) compared with the scatter of 

—TL correlation 



E. 



"y,iso,52 



= (13.8 ±2.0) ( 
0.0173 s 



E n 



513 keV 

0.67±0.43 



2.18±0.53 



(5) 



with a reduced x 2 = 2.28 (20 dof). The reason for the difference 
with Eq.|4]lies in the errors on T a , which are particularly large (see 
Tab. 1). To prove it, we have reduced the errors on T a by a factor 
1 .5 and recovered the result of the fit with a simple linear regression 
analysis (Eq.@. 

In Fig. [4] we compare the scatter of the data points around 
the E p — -E 7 ,i so correlation (defined by the 23 GRBs of Tab. 1) 
with the scatter of the data points around the best fit plane of Eq. [5] 
(i.e. obtained by weigthing for the errors). We note that the scatter 
of the data points in the i5 p -i? 7i i so -r a plane is comparable to the 
scatter of the data points in the E p -E lt i so plane. This suggests that, 
although in Eq.|5] (found by weighting for the errors) the exponent 
of T a is different from zero, the scatter of the data points is not 
improved with respect to the E p — E lt i so correlation. We conclude 
that it is Eq.|4] which better represents the best fit to the data. 

This would be sufficient to demonstrate that T a does not help 



E-y i so correlation. 



to define a tighter correlation between the isotropic energy and E p . 
To further demonstrate this point, we directly study the relation be- 
tween T a and E-,^ S o, comparing it with the corresponding relation 
between tj e t and i? 7i i so . 

Consider bursts with the same E p . For these the existence 
of the empirical Liang & Zhang correlation implies that _E 7i i so oc 
ij et -1 . This is proved by the data: in Fig. [5] we show the anticor- 
relation between tjet (rest frame) and E lt i so for the sample of 23 
GRB by considering bursts with similar E p (corresponding to the 
different symbols in Fig. [5}. The _E 7 ,i so oc £j ct _1 relation converts 
in i? 7 ,i so 6 l J 2 =t =const (the "Frail correlation") when using the relation 
between tj e t and 6>j (Eq. 1), for both the HM and the WM cases. 

Consider now Fig. [6] showing T a versus E lt i so : in this case 
T a does not anticorrelate with E lt i so , for constant E p . The absence 
of the anticorrelation translates in the impossibility to obtain a sim- 
ilar "Frail correlation" when using T a as a jet angle indicator (i.e. 
inserting T a instead of tj e t in Eq. 1). 
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and £break we find (under the assumption of small angles): 
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Figure S. The time tj et as a function of £ 7 ,i so for the burst of the G07 
sample. Bursts with different peak energy ranges (as labelled) are shown 
with different symbols. The dotted lines indicate tj ot cx: E~ iso , with differ- 

do follow this 



ent normalizations, 
law. 



to illustrate that bursts with the same E, 
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Figure 6. The time T a vs -E 7 ,i S o for the burst of this paper (Tab. 1). As in 
Fig.[5]we use different symbols for bursts with E p in different ranges. 



3.3 The origin of the E p — _E 7ia correlation 

Although T a spans 3 orders of magnitudes, and although it is not 
a jet angle indicator, the E p — E lt a correlation exists and has 
a scatter comparable (or slightly larger) than the scatter of the 
E p — -E 7 ,j so correlation. What is, then, its origin? 

In Eq. Q]the jet angle depends on both the break time (either 
T a or tjot) but also on i5 7 ,i so /»7 7 . This term is the kinetic energy of 
the GRB outflow left-over after the prompt radiative phase. Since, 
usually, one assumes a constant value 7/ 7 = 0.2, the collimation 
correction depends on E 1: i so . 

In this section we explore the role of E lt i so and tbrcak (either 
identified as T a or tjet) in Eq.^in shaping the Ep — E^ tdb and the 
correlation and their scatter. 

From Eq. Q] considering only the dependencies from -E 7 ,i so 



WM : £ 7 



-^7, iso OC E^y^is 



1/2 



HM : E-, oc E-, iso 9 Z oc S 7 ; 



' > ° E i/4 

-y,iso 



7 l/2 
J 7,iso 

3/4 

break xE - 



1/2 .1/2 
7, i so break 



.3/4 ^3/4 
7, iso break 



(6) 



Consider those GRBs that have the same E p , but different 
-E 7i i ao . Even if they have the same tbrcak, the use of Eq. [6] im- 
plies a clustering of the corresponding collimation corrected E- 1 : 
GRB with larger -E 7 ,i S o would have smaller 9, and therefore have a 
larger correction for their collimation (and vice-versa). 

We calculate 9 and correct for the collimation factor assuming 
that the break time (either T a or t- Jc t) is constant and equal for all 
bursts. The use of Eq. [6] with a fixed tbrcak gives E-y oc Ey? a 

in the WM case, and E-y oc E^ so in the HM case. We plot E p 
versus this estimate of E-, in Fig. |7j (for both the WM and HM 
cases). The empty star symbols define a correlation which is less 
scattered than the E p — £ 7i i so correlation (empty circles) in the 
HM and WM case (left and right panels in Fig.|7j. Moreover, the 
scatter reduction is somewhat more significant in the WM case with 
respect to the HM case. 

Note that this result has no physical meaning and it is only 

1/2 

due to the fact that we plot the points in the E p -E^ iso plane for 

the WM case (or E p -E^* o for the HM case) instead of in the 
E p -E lt i so plane. Up to this point the reduction of the scatter of 
the E p — i5 7 ,iso correlation, and the difference in the slope with 
respect to the Amati relation, is due to the dependence of Eq. [6] 
(used to derive the jet opening angle) on E lt i so . 

We can now see if the use of the real break time (either T a 
or tjet - upper panels and lower panels in Fig. [7] respectively) im- 
proves or not the found correlations. The results in the case of T a 
are shown in the upper panels of Fig.|7]for the HM and WM case. 
Although a correlation between E p and i? 7 , a appears (filled cir- 
cles), we find that the effect of T a is to increase the scatter of the 
data points with respect to the same correlation defined with T a 



constant. Note that the scatter of the found E v 



correlation 



and WM case). This demonstrates that the variable T a , if treated 
as a jet break time to estimate the angle and derive the collimation 
corrected energy, leaves unaltered (or even increases) the scatter of 
the E p — E lt i so correlation. This effect is clearly shown in Fig. [7] 
by bursts with typical rest frame peak energy in the range 100-1000 
keV which are the majority in our sample. 

In the bottom panels of Fig. |7j we show the results obtained 
with tjot. As discussed in Sec. 2, the sample of bursts used to define 
the Ep — _E 7 correlation only marginally overlaps to that used to 
define the E p — E~ u ^ correlation: this is the reason why also the 
Ep — .E 7 ,iso correlation in the lower and upper panels of Fig. |7J 
are different. If first we set ij ot =constant and equal for all bursts 
we have, as before, a reduction (empty stars) of the scatter of the 
Ep — .E 7 ,iao correlation (empty circles). But now, if we assign to 
each GRB in our sample its tj at found from the optical light curve, 
we see that the scatter of the data points is further reduced (note in 
particular this effect in peak energy range 500-1000 keV). 

From these results we conclude that: 

• the dependence of the jet angle on E lt i so reduces the scatter 
of the Ep — E lt i so correlation; 

• the use of the time T a which marks the end of the plateau 
phase in the X-ray light curves, increases the scatter of the correla- 
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Figure 7. Comparison of the scatter of the E p — S 7 , a and E p — E 1 correlations in the homogeneous (left) and wind (right) medium cases. Top panels: from 
right to left we report the E p — i? 7 ,i so correlation (empty circles), the correlation defined by assuming a constant value of T a (for all the bursts) in computing 
the collimation correction (empty stars) and the correlation defined by using the measured T a (filled circles; values of T a reported in Tab. 1). Bottom panels: 
same as above but considering tj ct instead of T a . 



tion found assuming a constant T a . This effect, counterbalanced by 
the previous one, makes the scatter of the E p — -E 7 , a correlation 
similar to the scatter of the E p — i5 7 ,i so correlation; 

• the (albeit weak) positive correlation between T a and -E 7l i so 
(see Fig.[6]l, makes the slope of the E p — -E 7 , a correlation to differ 
somewhat from the slope of the correlation found with T a constant; 

• the small dispersion of the Ghiiianda relation is due not only 
to the presence of -E 7l i so in the equation for the angle, but also to 
the jet break time ij e t . 



4 SUMMARY AND DISCUSSION 

We have analyzed in detail the correlation found by W07 between 
E p and the collimation corrected energy i? 7 , a , found using T a as 
a jet break. Our sample (of 23 objects) is larger than what used 
by W07, because we could add some GRB occurred after August 
2006, and some other bursts for which we could find the value of 
E p . We find that E p and -E 7 , a correlate, by we also find that this 



correlation is weaker, and has a different slope, than the E p — E~, 
one. 

Investigating further, we have also shown that the E p — -E 7 , a 
correlation is entirely a consequence of the existence of the 
Ep — i? 7! iso (Amati) correlation, coupled with the fact that the 
opening angle 8j depends on E lt i so . This is why one finds a corre- 
lation which is narrower than the Amati one even using the same 
T a for all bursts: in this case the found correlation is even better 
than the one obtained introducing the real T a . This demonstrates 
that T a does not play any role in the E p — -E 7 , a correlation: it is 
even worsening it. 

We have then made the very same test to the E p — _E 7 
(Ghirlanda) correlation, and found that in this case the use of 
the "real" £j et improves the correlation, decreases the scatter and 
slightly changes the slope of the correlation, with respect to the one 
obtained using a fixed ij e t- 

All these results are fully confirmed by analyzing the E p - 
£ 7 ,iso-T a correlation, and comparing it to the Liang & Zhang cor- 
relation (which uses ij ct instead of T a ). Again, all evidence is to- 
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Figure 8. Top panel: the time T' & as a function of the isotropic energy of the 
prompt emission E^ - mo calculated in the 15-150 keV energy range, with- 
out including any bolometric and fe-correction for all bursts with measured 
redshift (circles) and for those which also have a measured E p (stars). The 
insert shows the bolometric fluence vs the 15-150 keV fluence, to show that 
the two are proportional. Bursts in the insert are the ones in Tab. 1 . Bottom 
panel: the time T a as a function of the energy emitted in the plateau phase, 
in the 0.3-10 keV energy range. The fluences of the plateau phase have been 
taken from W07. Symbols as above. 



wards T a playing no role: the found E p -E lt i so -T^ correlation has 
a similar (if not larger) scatter than the E p — E lt i so correlation. 

We conclude that the time T a at which the plateau of the X- 
ray light curve ends is not an important parameter for defining the 
spectral-energy correlations. As a consequence, it is not a jet angle 
indicator. Instead, we confirm that tjet is likely to measure the jet 
opening angle, since it passed our simple test to see if it helps or 
not to strengthen the correlation between E p and the collimation 
corrected energy derived using a fixed value of tjet for all bursts. 

We note a weak correlation between T a and i5 7 ,i so when con- 
sidering the 23 GRBs in our sample (see Fig. [6] probability of 
chance correlation equal to 4 x 10 -3 , according to the Spearman 
rank correlation), which however disappears when considering the 
entire sample of GRBs of known redshifts and T a (and relaxing 
the requirement to have also E p known, see Fig. [8]>. Therefore we 
can also conclude that T a is not related to E lt i so , differently from 
tj e t which is instead anti-correlated with E lt i so when considering 
GRBs of similar E p . Note that, for a distribution of T a not cor- 
related (or anti-correlated) with i? 7 ,i ao , the resulting E p — £? 7>a 
correlation should become parallel to the one obtained fixing T a to 
a constant value. We therefore predict that, for larger samples of 
GRBs with known redshift, T a and E p , the E p — -E 7 , a correlation 
will have a slope « 0.75 (w 1) in the HM (WM) case. 

Finally, there is no correlation between T a and the isotropic 



energy (in the observed 0.3-10 keV energy range) emitted during 
the plateau phase (see Fig. [8}. 

Unfortunately, our finding that the E p — £? 7>a correlation is 
only a by-product of the E p — £/ 7j js (Amati) correlation is not 
very helpful in shedding new light on the problem of explaining 
why the light curve of most GRBs is characterized by the steep- 
flat-steep behaviour, and also why this is sometimes different by 
the simultaneous behavior observed in the optical band. 

Recently, we (Ghisellini et al. 2007) have proposed that the 
X-ray light curve is produced by internal dissipation in late shells, 
where "late" means that the central engine continues to create rela- 
tivistic shells even days after the trigger. This implies that the cen- 
tral engine is long lived, and that there is a sharp transition between 
the early phase (determining the duration of the early prompt, i.e. 
the time Tgrj), and the late phase, emitting much smaller luminos- 
ity, but for a longer time. The early phase should be characterized 
by a large bulk Lorentz factor Y of the created shells, changing 
erratically, while the late phase is instead characterized by shells 
created with Y decreasing monotonically. This implies that when 
r > l/f?j , we see only a fraction of the emitting surface, but when 
r ^ l/#j we see the entire emitting surface. This means that there 
should be a break in the light curve when Y — l/9j. We associate 
this time to T a . 

According to this idea, the two times T a and tjet have an ob- 
vious link: both are determined by the jet opening angle 8j . Then, 
why we do not observe a relation between these two times? We 
propose the following answer: the time tjet is determined when the 
fireball made by the merged shells of the early phase, carrying a 
kinetic energy E^, is decelerated by the external medium. The dy- 
namics of the process is robust, depending only on the value of the 
external medium density (and its profile with distance), on Ek, and 
finally on the conservation of energy and momentum. We have the 
same Y(t) law for different GRBs, if the external medium has the 
same profile. 

On the contrary, the law characterizing how Y changes as a 
function of time during the late prompt phase can be different from 
burst to burst. As far as we know, there might not be a robust analog 
of the conservation of energy and momentum dictating the Y(t) 
behaviour in this case. Indeed, we can have information on this 
law from the observed slopes of the light curve (see Ghisellini et 
al. 2007), and under simplifying hypotheses we can already derive 
that Y(t) is different from one bursts to another. As a consequence, 
there is no general formula, equal to all bursts, in which T a can be 
inserted to derive the jet opening angle. 
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